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1. Introduction 



By now there exist four well established ways of systematically constructing general 
solutions of the Einstein equations, by solving 

(i) a spacelike Cauchy problem (see [12, 1] and references therein), or 
i^l (ii) a boundary- initial value problem [24, 31, 23] (for further references see [3(S]), or 

(iii) a characteristic Cauchy problem on two transverse hypersurfaces, or 
\/J ■ (iv) a characteristic Cauchy problem on the light-cone. 

^^ ' One can further consider mixtures of the above. The aim of this paper is to present 

^^ , some new approaches to the last two questions, and to review the existing ones. 

To put things in perspective, recall that Einstein's equations by themselves do not 

have any type that lends itself directly into a known mathematical framework which 

\_j ' would provide existence and/or uniqueness of solutions [2"i]. The monumental discovery 

Z^ I of Yvonne Choquct-Bruhat in 1952 [i:»] was, that the imposition of wave-equations 

on the coordinate functions led to a system where both existence and uniqueness 

could be proved. The constraint equations satisfied by the initial data on a spacelike 

hypersurfacc turned out to be both necessary and sufficient conditions for solving the 

j^ , problem. The constraint equations, and the "harmonicity conditions" became then the 

JH ' two standard notions in our understanding of the spacelike Cauchy problem. 

- - -' In the early 1960's there arose a strong interest in the characteristic initial value 

problem because of attempts to formulate non-approximate notions of gravitational 

radiation in the non-linear theory [33, 2, 34, 37, 17]. While those papers provided much 

insight into the problem at hand, it is widely recognized that the first mathematically 

satisfactory treatment of the Cauchy problem on two intersecting null hypersurfaces is 

due to Rendall [3{)], see also [(>, 5, 14, 35, 21, IG, 32, 3, 4, 20, 27]. Kendall's initial data 

consist of a conformal class of a family of two-dimensional Riemannian metrics 

J '■— jAsirTX )dx dx 

on the null hypersurfaces, complemented by suitable data on the intersection. Here r 
is an affine parameter on the null geodesies threading the initial data surfaces. Rendall 
uses the Raychaudhuri equation to compute the conformal factor fl needed to determine 
the family of physically relevant data 

9 ■■= 9AB{r, x'^)dx'^dx^ = n'^{r, x'^)-fAB{r, x'^)dx'^dx^ (1) 

on the null hypersurfaces. The harmonicity conditions and the Einstein equations 
determine then uniquely the whole metric g to the future of the initial data surfaces 
and near the intersection surface S. The reader will find more details in Section 2. 
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Rcndall's elegant approach works well in vacuum, and more generally for a class of 
matter fields that includes scalar, Maxwell, or Yang-Mills fields. However, it appears 
awkward to use an unphysical family of conformal metrics as initial data, since the 
physically relevant, and geometrically natural, object is the family g. In this context it 
appears appropriate to view the tensor field g as an initial datum on the characteristic 
surfaces, with the Raychaudhuri equation playing the role of a constraint equation. 
The idea of prescribing 7 should then be viewed as a conformal ansatz for constructing 
solutions of this constraint equation. 

The last point is only a question of interpretation. More importantly, Rcndall's 
scheme does not work for e.g. the Einstein- Vlasov equations for particles with 
prescribed rest mass m (see [7] for an existence theorem for those equations with 
initial data on a spacclike hypersurfacc) , because the energy- momentum tensor for the 
Vlasov field,! 



Tc^p = Stt / fp^pp dfi{p) , (2) 

where / = J{x,p) is the Vlasov distribution function and d/i = d/^(p) is the Riemannian 
measure induced on the "mass-shell" {gpap''p'^ = —m^}, depends explicitly upon all 
components of the metric. This leads to the need of reformulating the problem so 
that the whole metric tensor is allowed as part of the initial data on the characteristic 
surfaces. Such a method will be presented in Section 3, after having reviewed Rcndall's 
approach in Section 2. We will do this both for data given on two transversally 
intersecting null hypersurfaces, and on a light-cone. 

We complement the above with a geometric formulation of the characteristic initial 
data in Section 4, where we give geometric interpretations of n, out of n-|- I, wave-map 
gauge constraint equations. 

The bottom line of our analysis is, that the gravitational characteristic initial 
data have to satisfy one single constraint equation, the Raychaudhuri equation. This 
raises the question, how to construct solutions thereof. In Section 5 we present several 
methods to do this. In the short Sections 5.1-5.4 we recall how this has already been 
done in the preceding sections. In Section 5.5 we analyse the Hayward gauge condition 
K = T/{n — 1), which may be used as alternative to an affine- parameterization-gauge 
where the function n vanishes. 

It has been proposed to use the shear tensor a as the free initial data for the 
gravitational field rather than 7. However, one defect is that it is not clear how to 
guarantee tracelessness of a. We present in Section 5.6 a tetrad formulation of the 
problem to get rid of this grievance. Finally, we adapt in Section 5.7 to any dimensions 
an approach of Helmut Friedrich (originally developed in dimension four, using spinors), 
where certain components of the Weyl tensor are used as unconstrained initial data for 
the gravitational field. Again, this requires to work in a null-frame formalism to take 
care of the tracelessness of the Weyl tensor. 

In the case of a light-cone the Hayward gauge leads us to the issue: Under which 
conditions is the assumption, that the vertex is located at the origin r = of the 
adapted coordinate system, consistent with regularity at the vertex? This question is 
considered in an appendix. 

2. Kendall's approach 

In this section we review Rcndall's approach to the characteristic initial value problem. 
For definiteness in the remainder of this section we will consider the vacuum Einstein 
equations; we comment at the end of this section on those energy-momentum tensors 
which are compatible with the analysis here. 

X We have included a factor Stt in the definition of Tij,^ , so that the Einstein equations read Si_iu = Tf^u, 
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Consider two smooth hypcrsurfaces Na^ a = 1, 2, in an (n+l)-diniensional manifold 
^, with transverse intersection along a smooth submanifold S. Near the N^s one can 
choose adapted coordinates {x^,x'^,x^) so that iVi coincides with the set {x^ ~ 0}, 
while N2 is given by {x^ = 0}. The hypcrsurfaces Na are supposed to be characteristic, 
which is equivalent to the requirement that, in the coordinates above, on A^i the metric 
takes the form 

g\N, = 5ii(dx^)^ + 23i2da;Mx2 + 2gi^da;Ma;'^ + g^^dx^da;^ , (3) 

= ■■9 

similarly on N2. Here, and elsewhere, the terminology and notation of [11] is 
used, in particular an overbar over a quantity denotes restriction to the initial data 
surface iVi U A^2- Rendall assumes moreover that x^ is an afSne parameter along the 
curves {a;^ = 0,x^ = const"*}, and that x^ is an affine parameter along the curves 
{a;2 = 0,a;'* = const-*}. 



On A^i let 



r^V^525^5 (4) 



2 

be the divergence scalar, and let 

aAB = ^d2'9AB - ';^~\^^AB (5) 

be the trace- free part of 925abi ^'Iso known as the shear tensor. The vacuum 
Raychaudhuri equation, 

d2T + \a\^ + ^=0, (6) 

n — 1 

provides a constraint equation on the family of two-dimensional metrics x^ i— > 
gj^g{x^,x'~^)dx'^dx^ , where 

\a\ ^ aA (Tb , cta = g ctac ■ 

Note that aA^ depends only on the conformal class of 'gAB- ^^ shown by Rendall, 
metrics satisfying the constraint (6) can be constructed by freely prescribing the family 
x^ I-)- ^AB(^x^Ax^ . Writing g^^ = 17^7^3, (6) becomes then a second order ODE in 
x^ for fi. 



Q = dln 



Vt 



ia2(7^^927As) + kP + ^^^[^^^d^lAB? 

2 4(n — 1) 



+-^l^''d2iABd2n ■ (7) 

n — 1 

This needs to be complemented by J7|s and 92^|s- 

Let us require all coordinate functions to satisfy the scalar wave equation DgX^ = 
0. Then the afhne parameterization condition rlgJAfi = can be rewritten as 

929x2 = -jT'9i2 ■ 

This equation determines the metric function 1712 on iVi with the freedom to prescribe 
512 on S. _ 

The equation R2A = on A^i takes the form 

- ^^2 + t)U + Vb^t^^ - ^a^r = , (8) 

where V is the covariant derivative operator of the metric g^sda; dx'^. Here, using 
the assumption that all coordinate functions satisfy the wave equation, the covector 
^A reads ([11, Equation (8.25)]) 

U := - '2g^^d2gxA + ^t^ViB^A"" + 2t'9iAr - 5As5''''f g^ , (9) 
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where the F^^i's arc the Christoffel symbols of the metric gj^gdx^dx^ . This provides 
an ODE for the metric functions gu^'. Indeed, one determines ^a by integrating (26), 
with the freedom to prescribe 

^^^ := U{^.' = 0) 

on S. (One should keep in mind that ^a here is unrelated to the corresponding field 
(,A on N2, determined by an analogous equation where all quantities r, a, etc., are 
calculated using the fields on N2.) 

Then g^A i^ found by integrating (9). Since the metric needs to be continuous, 
the metric component giA has to vanish at S; this requirement defines the integration 
constant. We further observe that by definition of (,a the freedom to prescribe ^^^ 
corresponds to the freedom of prescribing d2giA on S. 

The equation 'g Rab = in vacuum takes the form 

{d2 + r)c + i? - \g''''UiB + g^^'^A^B = , (10) 

where R is the curvature scalar of g, and where 

C:={2d2+T)f\ (11) 

Taken together, those equations provide a second order ODE for g^^; integration 
requires the knowledge of 'g^^ and <92ff^^ on S. Employing the relation g^^ = 
{'S^'^yiV^^ViAViB ^ 5ii) ■^e observe that g^^ has to vanish at S, while there remains 
the freedom of prescribing d2g^^, equivalently 925ii, on S. 

However, the validity of the harmonicity conditions implies certain constraints on 
S (see below): The value of d2gii at S is determined by equation (12c); similarly the 
function 92ffii at S follows from (12a). 

One has thus determined all the metric functions g^^ on iVi; the procedure on 
N2 is completely analogous. These are the data needed for the harmonically-reduced 
Einstein equations, which form a well-posed evolutionary system for the metric. 

However, not every solution of the equations constructed in this way will satisfy 
the vacuum Einstein equations: One still needs to make sure that the harmonicity 
conditions are satisfied. There is in fact one more subtlety, as one needs to verify that 
the parameter x^ is indeed an affine parameter on the null geodesies threading A^i. It 
turns out [36] that all this will be verified provided three more conditions are imposed 
at S: If we write v^ for what was giA\Ni so far, v^ for 32^1^21 the wave-coordinates 
conditions will hold if we require that on S 



p. I _ d2yJdeXgAB .-,„ . 

Vdct gAB 

1 



div^ + ^2i^A = /J . 9ABdc{gi2\/'^ctgEFg^^) , (12&) 

Vdet gEF 

p. I _ di^ del gAB ,,„ , 

Vdct gAB 

As already indicated above, the integration functions dig22\s and 925ii|s cannot 
be freely specified but have to be chosen in such a way that the equations (12 a) and 
(12c) are fulfilled. Equation (12&) will be satisfied by exploiting the freedom in the 
choice of div~^ and d2v\. So there remains the freedom to prescribe, say, div^ — d2v\. 

The constraint equation (126) can be tied to a terminology introduced by 
Christodoulou [1.3] as follows: Let L and L be two null normals to a codimension- 
two spacelike hypersurface S satisfying 

5(L,L) = -2. 

Christodoulou [13] defines the torsion one-form of S by the formula 

ax) = \g{VxL,L) . (13) 
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where X <E TS. Assuming gi2|s is positive we can choose L ~ \/2g^^ 82 then, on S, 
using the notation above, L = —\/2§^di^ and (13) reads 

= \t\di92A - d29^A) = Iv^'idii^A - d2^A) ■ (14) 

So C,A contains precisely the information needed to determine dii^^ and d2i']j at S, 
after taking into account {12 b). 

Theorem 2.1 (Rendall) Consider two smooth hypersurfaces Ni and N2 in an (n + 
1) -dimensional manifold with transverse intersection along a smooth submanifold S 
in adapted null coordinates. Let ^ab be a smooth family of Riemannian metrics on 
Ni U A'^2; continuous at S . Moreover, let 57, difl, 82^, f and fA, A = 3, . . . , n + 1, 
be smooth fields on S , where we assume n and f to be nowhere vanishing on S. Then 
there exists an open neighbourhood U of S in the region {x^ > 0, a;^ > 0}, a unique 
function D, on (Ni U N2) U and a unique smooth Lorentz metric g^^ on U such that 

(i) g^i, satisfies the vacuum Einstein equations, 

(^i) 9ab = ^'^lAB, 
(Hi) fl induces the given data on S, gi2\s — f o,nd C,A ^ fA- 

This analysis of the constraints applies equally well to a light-cone with some 
minor modifications [11], where the wave-equations for the coordinate functions are 
replaced by wave-map conditions. Moreover, there is no need to provide further initial 
data at the tip of the light-cone, as those are replaced by conditions arising from 
the requirement of regularity of the metric there; the reader is referred to [11] for a 
detailed discussion. An explicit parameterization of tensors g which arise by restriction 
of a smooth metric in normal coordinates has been given in [15]. The reader should 
keep in mind the serious difficulties with regularity of the metric at the vertex, when 
attempting to prove an existence theorem for the light-cone problem; see [10, 9] for 
results under restrictive conditions on the data. 

The above extends easily to non- vacuum models with energy-momentum tensors 
of the form 

T22 = T22(matter data,^AB,dijAB,^,d2fl,gi2,digi2,x^) , « = 2, A , 
T2A = T2A(matter data,7AB,9i7AB,f^,9jf2,5i,,9i5;^2;^2ffiAit^iff22,a;') , 
T12 = Ti2(matter data, 7AB,<9j7AB,i^, SiSl, 31^,9231^, 9^312, c?ig2j, a;') , 
on the initial surface {x^ = 0}, cf. [11]. 

3. All components of the metric as initial data 

Let i'^ denote the field of null tangents to a characteristic hypersurface. In this section 
we present a treatment of the characteristic Cauchy problem which applies to energy- 
momentum tensors of the form 

r^X = Tp(.9,0,all,g,all0,a;), (15) 

for some fields (f> satisfying equations which, when the metric is considered as given, 
possess a well-posed characteristic Cauchy problem. Here the symbol c?" denotes 
derivatives in directions purely tangential to the initial data surfaces. In particular 
(15) includes the Einstein- Vlasov case. 

As already discussed, in [-Mi] the corresponding problem for the vacuum Einstein 
equations is solved using an affine parameterisation of the generators and a wave- 
map ("harmonic") gauge. In Kendall's approach some components of the metric are 
calculated by solving the characteristic-harmonic gauge constraint equations, which 
form a hierarchical ODE-system along the generators of the initial surface. For an 
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energy-momentum tensor (2) this approach will generally lead instead to a quasi-linear 
PDE-system for the metric components. To establish an existence result for that system 
might be intricate, if possible at all. It is in any case not obvious how to include an 
energy- momentum tensor (2) in this scheme, compare [.39]. We circumvent the problem 
by using a gauge adapted to the initial data, where the metric tensor, and thereby (2), 
is fully given on the initial surface, while the wave-gauge source vector W^ is computed 
from the values of the metric on the initial surface using the Einstein ~ wave-map- gauge 
constraint equations of [11]. 

We start with an analysis of two intersecting hypersurfaces, the case of a light-cone 
will be covered in Section 3.2. 

3.1. Two transverse hypersurfaces 

Consider two smooth hypersurfaces Na, a = 1, 2, in an (n -I- l)-dimensional manifold 
^, with transverse intersection along a smooth submanifold S. As before, we choose 
adapted null coordinates (x^,x^,a;"*) so that A^i coincides with the set {x^ = 0}, while 
N2 is given by {x^ = 0}. We use a "generalized wave-map gauge" as in [11], with 
target metric g of the form 

g = 2da;M2:^ -^ gAB{x^,x'^, x^) dx^dx^ . 

Here gAB is any family of Riemannian metrics on S parameterized by x^ and x^, 
smooth in all variables. The metric g is only introduced so that the harmonicity vector 
iif, defined in equation (18), is a vector field, and plays no significant role in what 
follows. 

As gravitational initial data on the initial hypersurfaces we prescribe all metric 
components g„y in the coordinates above, as well as a connection coefficient k; this 
needs to be supplemented by initial data 4> for (p. For instance, in the Einstein- 
Vlasov case, the supplementary data will be a function / defined on the mass-shell 
{'SiivP^P'^ = — m^}, viewed as a subset of the pull-back of T^ to the A^q's. 

The hypersurfaces Na are supposed to be characteristic, which is equivalent to the 
requirement that, in the coordinates above, on A^i the metric takes the form 

9\N^ = 9iMx^f + 2gi2da;Ma;2 + 2gi^da;Mx'^ 4- 5ABda;'^dx^ , (16) 

=■■9 
similarly on N2. Here, and elsewhere, the terminology and notation of [11] is used, 
in particular an overbar over a quantity denotes restriction to the initial data surface 
Ni U A^2- (Some obvious renamings need to be applied to the equations in [11], for 
instance the variable u there is x^ on A^i, and x^ on N2] the variable r there is x'^ on 
Ni and x^ on A''2.) 

We want to apply Rendall's existence theorem [.■^t(]] for an appropriately reduced 
system of equations. For this the trace, g^^,, of the metric on the initial surface A'^i U A^2 
needs to be the restriction of a smooth Lorentzian spacetime metric. This will be the 
case if gj^j is nowhere vanishing, if 'gAB\Na is ^ family of Riemannian metrics, and if 
g„j^ is smooth on A^i and A'2 and continuous across S = Ni O N2. We therefore need 
to impose the following continuity conditions on S, 

lim gAB\Ni = lim .gAsJTVa , (17a) 

lim gi2\Ni = lim .gialA^a , (17&) 

lim giA\Ni = , lim g2A\N2 = , (17c) 

lim giilA^i = , lim 522|A^2 = . (17d) 

x2— j-o aji-i-O 



Let iJ^ be the harmonicity vector, defined as 

^ap - W^ , with W'.^y X „^ 



H^ := g'^'^T^p - W^ , with W^ := 5"^f ^. +W^ , (18) 
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where W'^ will be a vector depending only upon the coordinates, and determined by 
the initial data in a way to be described below. (In principle W^ can be allowed to 
depend on the metric as well, but not on derivatives of the metric.) To obtain a well 
posed system of evolution equations we will impose the generalized wave-map gauge 
condition 

H^ = . 

More precisely, we view the wave- map gauge constraints [11] as equations for the 

restriction W^^ of W^^ to A^i and N2. We will solve those equations hierarchically. 
We emphasize that, assuming (15), all components of the energy momentum tensor 
restricted to A^i and A'^2 are explicitly known since g„y and (j) are. 

We present the calculations on iVi, the equations on N2 are obtained by 
interchanging the index 1 with the index 2 in all the formulae. 

Let 5*^11/ denote the Einstein tensor. In the notation and terminology of [11] the 
first constraint, arising from the equation S22 = R22 — T22 evaluated on iVi, reads 
(see [11, Equation 6.11]) 

^2 _ 

-d2T + KT-\a\^ -=T22, (19) 

n — 1 
where r and a are defined as in (4) and (5), respectively. Indeed, using the formulae 
in [11, Appendix A] one finds 



w^n, = d2Ti2 + {Ti2r - TI2TI2 



S22 ~ diTl^- d2(T^2-^2A) + (Xi2+^2a)^22 - (^12) ^^23^2^ 

— - O2I 2yl + ^ 2A^ 22 ~ ^ 25^ 2A 

^ -d2T + TTI2 - xa^'xb'' , (20) 



where 



XA^ := ^g^'^fts^c 



Here we adapt the point of view that the function k is the value on Ni of the 
Christoffel coefficient T22, and is part of the initial data. Hence, we view (19) as a 
constraint equation linking g^^, k, and the matter sources (if any). 

In the region where t has no zeros, (19) can be trivially solved for k to give 

T 

It follows that K does not need to be included as part of initial data when r has no 
zeros, and then the constraint equation (19) can be replaced by the last equation, 
determining k. 

Equation (21) can still be used, by continuity, to determine k on the closure of 
the set where r has no zeros, keeping in mind that the requirement of smoothness of 
the function so determined imposes non-trivial constraints on the right- hand- side. In 
any case (21) does not make sense if there are open regions where r vanishes. It seems 
therefore best to assume that k is any smooth function on TVi such that (19) holds, and 
view that last equation as a constraint equation relating k, 'g^g and its derivatives, 
and the matter fields; similarly on N2. 

It should be kept in mind that, once a candidate solution of the Einstein equations 
has been constructed, one needs to verify that k, is indeed the value of T22 on A^i . We 
will return to this in (39). 

We choose W^ to be 

W^ ■= -W^- t^{2K + r) - 2^25^2 (22) 

Note that the right-hand side is known, so this defines W^ . By definition, this is 
the di component of VF^ in the coordinate system {x^.x^.x^). We will define the 



The many ways of the characteristic Cauchy problem 8 

remaining components of W^ shortly, the resulting collection of functions transforming 
by definition as a vector when changing coordinates. 
From [11, Appendix A] one then finds 

rL = « - ImH' , (23) 

and (20) together with (19) give 

'S22-T22 = -lmH^r. (24) 

The corresponding constraint equation on 7V2 determines PF^lAfg. We shall return 
to the question of continuity at S of W^\n^dn,^ and of W^\i^^ijN2 shortly. 

The next constraint equation follows from S2A = R2A = T2A- From the formulae 
in [11, Appendix A] we find 

^1^2/1 = dA^i2 + ri2(ri^ — r2^) + ^12^ AB ^ ^2A^1B ' 

which gives 



S2A = diT\j^ + a2r2^ + dsT^A - dAV\2 - OatIj -dAris + r^sPf^^ 

I "pi /'"p2 T^l \ I p-B "p2 I "pi? j^C "pi j^B j^B "pO 

+^ 12l^ 2A~ '^ lAI + ^ 2B^ 2A + ^ BC^ 2A ~ ^ AB^ 12 ~ ^ AC^ IB 

= ^2r2^ + ^sFa^ — 9aF22 - dA^2B + r2_Br2^ + YbC^2A ~ ^AC^2B 

= {d2 + r)Tlj, + VeXa"" - dATl2 - Oat , (25) 

where V is the covariant derivative associated to the Riemannian metric qab- That 
leads us to the equation 

- ^(^2 + t)U + Vb^t^^ - "^dAT ~ Oa^ = T2A , (26) 

z n — 1 

where the field ^a denotes the restriction of the (rescaled) ChristofFel coefficient — 2F2^ 
to A'^i. We determine ^a by integrating (26), with the freedom to prescribe 

on 5. (One should keep in mind that ^a here is unrelated to the corresponding field 
(,A on JV2, determined by an analogous equation where all quantities r, a, etc., are 
calculated using the fields on N2) We then define W^ through the formula 

W^ := g-^^ [^B + 2t\d2g,B - ^Vic'^b'' - Vibt) - ViBiW' + W') 

+ g''''r^D^W^; (27) 

equivalently 

U = -29'^d2g^A + 45''5iB^A^ + 2t%AT + g.^W' + W') 
+9AB {W'' + W'')- gAsf'^cD ■ (28) 

This has been chosen so that, using the formulae in [11, Appendix A and Section 9], 

S2A^T2A = -\{d2+T){gABH'' +giAH') + \dA{-gi2H^) ■ (29) 
Moreover, one finds (cf. equation (10.35) in [11]) 

U - ~2TIa - VabH"" - ViaH' . (30) 

On S (27) takes the form 

717^1 —AbI^Ni . r,— 120 — 1 I —BCtrA pA N 

W \s = g I^Cb +25 025ibJ+5 (^bc~^bc)- 
Keeping in mind the corresponding equation on A^2, 

W \s = g ^B +2g oig^B +g (^bc-^bc)- 
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the requirement of continuity of W^ |Ariu7V2 leads to 

^a' - ^a' = 29'\d,g2A - d2giA)\s ^ <a ■ (31) 

Recall that the torsion one- form C,a has been defined in (14). 
We continue with the equation S12 = T12, or, equivalently, 

ff^^i?AB = -{'it^T^2 + t^T22 + 2f^T2A) = V'^^'Tab - T , 

which we handle in a manner similar to the previous equations. Using the identities 
(10.33) and (a corrected version of §) (10.36) in [11] we find that on TVi we have 

g^^RAB ^ {d2 + rL + r)(25^^ri5 + rf') ^-g^^Vl^Tl^ - 2g^^V^r^^ + R , 

and we are led to the equation 

{d2 + K + t)C + (V^ - \u)i^ + R = V^^'Tab - T , (32) 

with S^^ := g^^^B, and where the quantity C denotes the restriction of 

^[9 '-AB + T^g ) 

to Ni. We integrate (32), viewed as a first-order ODE for C. The initial data on S are 
determined by the requirement of continuity of W^ at S*, which we choose to be 

W^--^lc-{d2 + n+^T)f^-W^ . (33) 

Indeed, recall that W^ \ 5 has already been calculated algebraically when analysing the 

first constraint equation on N2, in exactly the same way as we calculated W^ in the 
first step of the analysis above. 

Similarly the initial data for the integration of the constraint which determines 
H^^Iat^ is determined by the requirement of continuity oiW^\NiUN2- 

The choice (33) has been done so that 

g^^'RAB - g^^'TAB +T = {d2 + K + r- ^g,Jl'){2W - ^,2t^W) - \cg12H' 

+ (Va -U- IvabH"" \g,AH'){H^ + g^cV^'^H') ■ (34) 

We note that our choice of W^ is equivalent to 

C = 25^^riB + rt' + g,2-f^W - 2H^ . (35) 

Summarising, given the fields k, g„j, and (j) on TVi U N2, satisfying (31), and the 

sum ^^^ + '^A^ o'^ "^i ^® have found a unique continuous vector field W on iVi U N2^ 
smooth up-to-boundary on Ni and N2, so that (24), (29) and (34) hold on iVi U A^2- 

Letting W be any smooth vector field on ^ which coincides with W' on iVi U 7V2, and 
assuming that the reduced Einstein equations (see (43) below) can be complemented 
by well-posed evolution equations for the matter fields, we obtain a metric, solution of 
the Cauchy problem for the reduced Einstein equations in a future neighbourhood of S. 

However, the metric so obtained will solve the full Einstein equations if and only 
if [11] H^^ vanishes on NiVJ N2; so we need to ensure that this condition holds. Note 
that at this stage a smooth metric g, satisfying the reduced Einstein equations, and a 
smooth vector field W^ are known to the future of Ni U N2 in some neighbourhood of 
5, and thus H^^ is a known smooth vector field there. 

By [11, Section 7.6], H^ will vanish on TVi if and only if H^ vanishes on S. Using 
(22) and (23) together with the equations in [11, Appendix A] we find 

H% = {g^^fd^g22 + 2ffi2^ + 2925I2 (3g) 

§ On the right-hand-side of (10.36) in [11], in the conventions and notations there, a term Tg^^/2 is 
missing. 



The many ways of the characteristic Cauchy problem 10 

We conclude that H^\n-^ will vanish if and only if the initial data ^32 '^^ ^2 have the 
property that the derivative dig22 on S satisfies 

dig22\s = 2(a2.gi2 - 912K.N1) <=> ^l2\s = i^Ni , (37) 

where, to avoid ambiguities, we denote by K7v„ the function k, associated with the 
hypersurface Na, etc. Similarly i^^|Ar2 will vanish if and only if we choose gn on A'^i 
so that 

92gii|s = 2(9i5i2 -gi2KAr3) ■^=> r}i|s = «;jv2- (38) 

With those choices we have H'^\s = H'^\s = 0, and the arguments in [1 i] show that 
H^\n2 ^ H^\N^ =0 as well. 

We continue with H . Then by equation (30) we have 



"^^(9.1512 - 9ig2A + d29iA) + 9abH^ + giAH^) 



C^^ = - [2TIj, + gABH'' + giAH^)\s 

e^^ = - {2T\a + gABH"" + giAH^)\s 
--{9 



^^'^{dAgi2 - d2giA + dig^A) + gAsH^ + giAH^) 



and the conditions H \s = and H^\s = determine ^^° in terms of the remaining 
data. Note that (31) is then automatically satisfied, and that we loose the freedom to 

prescribe Ca'+Ca '• 

It remains to show that our choice of the parameterization of the null rays is 

consistent, i.e. we have to make sure that the relations r22|Afi = uni and F^j^Jat^ = kjv2 

hold. This follows trivially from the vanishing of the wave gauge vector H due to the 

identities 

H^\N, = 2g^^{K^Tl^) and H^\N, = 2g^\K~T\^). (39) 

Similarly, the vanishing of H^ and H^ shows via (30) that the identification of S,a with 
the rescaled Christoffel coefficient — 2r2A on iVi and — 2rJ^ on N2 is consistent. The 
vanishing of H^ and H^ , together with the identity (35) imply that on A'^i the field 
C, indeed represents the value of 2'g^^T\g + rg^^, and the corresponding field on N2 
represents the value of 2g^^T\g + rg^^ there. 

In particular, the above provides a new and simple integration scheme for the 
vacuum Einstein equations, where all the metric functions are freely prescribable on 
iVi UiV2: 

Theorem 3.1 Given any continuous functions (^,3^1^) o?i A'^i U N2 such that 

5ki == 522(da;^)^ + 23i2da;Mx^ + 2g2A<^x'^<^x^ + SAsdx^da;^ , (40a) 
5^2 = 5ii(da;^)^ + 2:gi2<^x^<^x^ + 2giAda;Ma;^ +5A5da;^da;^ , (406) 

smooth up-to-boundary on Ni and N2, and satisfying (37)-(38) together with the 
vacuum constraint equations (here kn^ '■= i^\ni, etc.) 

T-2 

- d2TNi + KwiTATi - Ictatj P ^ = on TVi , (41a) 

n — 1 
2 

-Si Tat, + KN2TN0 - Wn2 P h- = on N2 , (416) 

n — 1 

there exists a smooth metric defined on some neighbourhood of S , solution of the vacuum 
Einstein equations to the future of Ni U N2 . 

Note that all the conditions are necessary. To see this, let g be any metric solving 
the Einstein equations to the future of iVi U 7V2, with Na characteristic. We can 
introduce adapted coordinates near A^i U N2 so that (40 a)- (40 6) hold. The constraints 
(41a)-(416) follow then from the Einstein equations [11], while (37)-(38) follow from 
our calculations above. 
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Proof: While the main elements of the proof have already been given, to avoid 
ambiguities we summarize the argument: Let {K,g) be given as above. Set 

^ := [0, oo) X [0, oo) X S* , 

where the first [0, oo) factor refers to the variable x^, and the second to x'^. On ^ let g 
be the metric g = 2dx^dx^ + (pAsdx dx^ , where (j)ABdx dx^ is a Riemannian metric 

on S. Let VV^ be constructed as above. Let W^ be any smooth extension of W^ to 
^, and let g be the solution of the wave-map reduced Einstein equations R^a = 0, 
with initial data y, where 

RiV ■= R^P - \{9c.xDpH^ + gpxb^H^), (42) 

with _ff^ defined by (18), and where D is the Levi-Civita covariant derivative in the 
metric g. (It follows from [s, page 163] that R)^o is a quasi-linear, quasi-diagonal 
operator on g, tensor-valued, depending on g, of the form 

Rif = -~\g^''D^b^g^p + f[g, DgU , (43) 

where f[g, Dg]ap is a tensor quadratic in Dg with coefficients depending upon g, g, W, 
DW and DW] existence of solutions of this problem follows from [Ui].) As H = by 
construction, a standard argument shows that H^ = 0, and so g is a solution of the 
vacuum Einstein equations in a suitable neighbourhood of S in ^. D 



3.2. The light- cone 

Now let us consider the same problem on a light-cone Co with vertex O. We prescribe 
the metric functions g„j^ on the cone in adapted null coordinates (cf. [11]) as well as (j) 
and, if r has zeros, k (note that t = \'g dig ^^ has no zeros in a sufficiently small 
neighbourhood around the vertex). In this section the notations and conventions from 
[11] are used again; in particular the x^-coordinate will be frequently denoted by r, 
and the light-cone is given as the surface {a;° = m = 0}. 

For Co to be a characteristic cone we need to have g^i = = Hi a ii^ our adapted 
coordinates. To end up with a Lorentzian metric, the component i^o = 5oi ^'"^^ ^'^ ^^ 
nowhere vanishing, while 'gAB has to be a family of Riemannian metrics on S*"^^. We 
consider initial data which satisfy 

ffoo = -1 + 0{r^) , 

z/o = l + 0(r2), 

UA = 0{r^) , 

9ab = r^SAB + 02{r^) , 

dldcdo'gAB = ^dcdoSAB 

for small r, where / = Onir") means that d^d'^f 
Sab denotes the round sphere metric. 

These conditions ensure that the metric is of the same form near the vertex as 
in [11]. The assumptions concerning the derivatives, which are compatible with the 

relations (4.41)-(4.51) in [11], are made to compute the behaviour of W near the 
vertex: 1 1 Though we do not attempt to tackle the regularity problem at the vertex 

here, as a necessary condition we want to make sure that W remains bounded near the 
vertex, which in our adapted coordinates means 



drVoo 


= 0{r) , 






(44 a) 


drVo = 


= 0{r) , 






(446) 


drVA -- 


= O(r^) 


, 




(44c) 


drdc'g. 


AB = 2r 


dcSAB + Oi{r^ 


) , (44d) 


+ 0{i'' 


'■), 






(44 e) 


^d^f- 


: O^r"^-^] 


)for 


i+|/3| <n. 


The tensor 



W°^0{1), W^^O{l), W^^O{r~^) 



II It is conceivable that a larger class of initial data turns out to be compatible with regularity at the 
vertex. 
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In fact it turns out that with (44 a)- (44 e) and the subsequent assumptions on the 

target metric and the energy momentum tensor the vector W goes to zero. 

We present the scheme for an arbitrary target metric g that satisfies the relations 

i>o = 1 + Oi(r2) , VA = Oi(r3) , §^„ = -1 + ©(r^) , (45a) 

dr9oo = 0{r) , Qab = r^^AB + Oi{r% (456) 

MTi = 0(r) , d^§^ = 0(r2) , g^'^MIE = 0{r) . (45c) 

Again, these assumptions are to ensure that the behaviour of W can be determined at 
the vertex. 

Additionally, we take a look at two particular target metrics: a Minkowski target 
5 = ry as in [1 1], and a target metric g ~ C which satisfies C ^g and which simplifies 

the expressions for the components of W. 

Let us now solve the constraint equations. The first constraint yields (supposing 
that T has no zeros, the case where it does have zeros can be treated as in the case of 
two transversally intersecting hypersurfaces) 

. . *i±iEi:!±tP±I;i (46) 

T 

and 

W°= -W°- i/°(2k + r) - 29iiy° . 

If we assume ^ 

Tn = 0(1) , dATn = 0(1) , dAdsTii = 0(1) , 

we obtain with our assumptions (44 a)- (44 e) and with the assumptions (45 a)- (45 6) 
concerning the target metric 

K = 0{r) , Oak = 0{r) , dAdsK = O^r) , 

and 

W° = 0{r) . 

Let us write = for an equality which holds when g is the Minkowski metric, with 

c 
an obvious similar meaning for =. Then 

#0 1 -rr^'sAB , 

and also 

M^o ^ 2i^"(f}, - k) . 

From the second constraint equation one first determines S,a- Recall that this 
is a first-order ODE. The integration constant which arises is determined by the 
requirement of finitcness of S,a at the vertex (cf. [1 1 , Section 9.2]), 

^bcja'' - "^—^Qat - Oak - Tia] df . (47) 

V n — 1 J 

If we assume 

TiA = Oir) , ObTia = Oir) , 
and employ (44 a)- (45 c) we find 

U = Oi(r2) . 

^ These assumptions on the energy-momentum tensor, as well as those which will be made later, will 
hold for a tensor T^,^ which has bounded components in coordinates which are well behaved near the 
vertex; note that the [u,r,x^) coordinates are singular at the vertex. 
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The function W^ can then be computed algebraically, 

= 0(1), 



where 



In particular 



and 



B _ 1_SCq — 

XA = 2^ Oi9ac 



W^l^y-gAB-g^^^-gBCgA^^ 



The functions Sg^j denote the Christoffel coefficients associated to the round sphere 
metric. 

Finally, we have a first-order equation for 

C = (25i + 2k + T)g" + 2W^ + 2W^ . (48) 

It reads 

{di + n + r)C + i? + ^'^^(VaCb - \uiB)+t^Tii + 2gi^TiA + 2i."Toi = . (49) 
This can be integrated. 



c 



ipll 



r^~-ie/:(«+-^)d?(^ 



-g^^'VACB - \g^''UiB+t^Tii + 2t^TiA + 2i^°Toi )dr 



where c is an integration constant. 

Using again the relations (44 a)- (44 e) and our assumptions on the target metric we 
deduce that 

g" = l + 0(r2), d,t'=0{r), 

R = {n - l){n - 2)r-'^ + 0{l) . 
Assuming that Toi = 0(1) we find that a general solution Q has a term of order r~("~^^ 
due to which W^ would not converge at the vertex. We thus set c = 0. That yields 

C = -(n-l)r-i+0(l) . 
Inserting this result into (48) we end up with 

W'^ = 0{r) . 
For that we employed 

W^ = -{n-l)r-^ +0{r) . 
In the special case of a Minkowski target we have 

Moreover, we find 

w'^lc-g^'^r\B-lrt'+nrl,-.). 

Let us assume that the vector field W^ can be extended to a smooth spacetime 
vector field W^ on the space-time manifold ^. If we further assume, as in the case of 
two transvcrsally intersecting null hypersurfaces, that the reduced Einstein equations 
can be complemented by well-posed evolution equations for the matter fields, for 
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sufRcicntly well behaved initial data we obtain [18] a solution of the Cauchy problem 
in a future neighbourhood of the tip of the cone. The metric obtained this way solves 
the full Einstein equations if and only if H vanishes on Co, as shown in Sections 7.6, 
9.3 and 11.3 of [11]. 

Let us assume now that initial data (k,5„i^) and a target metric g have been 
specified. In order to prove that the wave-map gauge vector H^ vanishes on the cone, 
one first establishes that it is bounded near the vertex. In our adapted coordinates 
that means 

i?" = 0(l), i?i=0(l), H^^Oir-^). (50) 

If we assume that those transverse derivatives which appear in the generalized wave- 
map gauge condition H ^ satisfy (compare [I I] for a justification under the conditions 
there)"*" 

d^,=0{r), d^^=0{r^), ff^^a^^ = 0(r), 
and the initial data fulfill, additional to (44a)-(44e), the relations, 

dAiyo = 0{r^) , dB<^A = 0{r^) , 
then one finds (using (44 a)- (45 c)), say in vacuum, 

i?" = Oir) , H^ = Oir) , H^ == 0(1) , 
which more than suffices for (50). 

4. A geometric perspective 

4-.1. One constraint equation 

Let us present a more geometric description of initial data on a characteristic surface. 

A triple (^, g, n) will be called a characteristic initial data set if ^ is a smooth n- 
dimensional manifold, n > 3, equipped with a degenerate quadratic form g of signature 
(0, +, . . . , +), as well as a connection form k on the one-dimensional degeneracy bundle 
Kerg, understood as a bundle above its own integral curves. The data arc moreover 
required to satisfy a constraint equation, as follows: 

We can always locally introduce an adapted coordinate system where Kerg is 
Spani9i.* (There only remains the freedom of coordinate transformations of the 
form {x^,x"^) !->■ {x^ {x^ , x^) , x^ (x^)) .) Then the connection form k reduces to one 
connection coefficient: 

ya^^l^K^l. (51) 

In this coordinate system we have g = 'g^giix'^dx^ . Denoting by y^^ the matrix 
inverse to g^g, set 

XB^:=lg^''d,gcB, t:^xa^. (52) 

Under redefinitions of the adapted coordinates the field r transforms as a covector, 
which leads to the natural covariant derivative operator 

ViT := {di - k)t . (53) 

With those definitions, the characteristic constraint equation reads 

Vir = -xs^XA''-p ^^ Vir+^-T^^-as^a^^-p, (54) 

71—1 

+ Note that these transverse derivatives are obtained from the solution g of the reduced Einstein 
equations with initial data <?. The assumptions (50) are known to hold e.g. if one uses the wave-map 
gauge W^ = near the vertex [11]. 

* This vector was denoted as dr or 9i in Section 3.2, by 9i in Section 2 when considering the null 
hypersurface N2, and as 82 in Section 2 when considering the null hypersurface A'^i. 
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where p represents the component Tii\^ of the energy- momentum tensor of the 
associated space-tniie (^, g) and, as before, a is the trace-free part of x- 

A triple {^,g,K) satisfying (52) with p = wih be called vacuum characteristic 
initial data. 

An initial data set on a light-cone is a characteristic data set where ^ is a star- 
shaped neighbourhood of the origin in M" from which the origin has been removed, 
with the tangents to the half-rays from the origin lying in the kernel of g, and with 
{g,K) having specific behaviour at the origin as described e.g. in [11]. 

The reader is referred to [2G] for a clear discussion of the geometry of null 
hypersurfaces, and to [29, 28, 30] for a further analysis of the objects involved. 

4-2. Dim-jV constraint equations 

An alternative geometric point of view, closely related to that in [SU], is a slight 
variation of the above, as follows: Instead of considering a connection on the degeneracy 
bundle Ker^, viewed as a bundle over the integral curves of Ker^, one considers a 
connection on this bundle viewed as a bundle over jV . For this one needs the connection 
coefficients k and ^a-, defined by the equations: 

Va ai = K^i , Va^ ai = - i^Si + XA^'Ss . (55) 

The coefficient k satisfies the same constraint equation as before. The remaining 
coefficients i^A are obtained from (26), in notation adapted to the current setting: 

- k^i + T)iA + ^bOa"" ' "^^dAT -dAK = TiA- (56) 

2 n — 1 

The fact that this system of ODEs can be solved in a rather straightforward way 
(compare (47)) given k and gAB, should not prevent one to view this equation as a 
constraint on the initial data. 

A useful observation here is that in [.30, Appendix C], where it is shown that 
(54) and (56) can be obtained from the usual vector constraint equation on a spacelike 
hypersurface by a limiting process, when considering a family of spacelike hypersurfaces 
which become null in the limit. 

The alert reader will note that the constraint equations (54) and (56) exhaust all 
tangential components of T^^l'^ . We are not aware of a geometric interpretation of the 
equation (49), which involves the remaining, transverse, component of Tf^ui'^ ■ Pursuing 
the analogy with the spacelike Cauchy problem, one could be tempted to think of this 
equation as corresponding to the scalar spacelike constraint equation. However, this 
analogy is wrong since it is shown in [30, Appendix C] that the scalar constraint 
equation and one of the vector constraint equations degenerate to the same single 
equation when a family of spacelike hypersurfaces degenerates to a characteristic one. 

4-.3. Uniqueness of solutions 

Given a vacuum characteristic data set on a light-cone, or two vacuum characteristic 
data sets with a common boundary S (where some further data might have to be 
prescribed, as made clear in previous sections), one can impose various supplementary 
conditions to construct an associated space-time metric. For example, one can redefine 
x^ so that K = and impose wave-coordinate conditions, or wave-map coordinate 
conditions in the light-cone case, to obtain the required space-time metric. Or one can 
prescribe the remaining metric functions as in Section 3, with appropriate conditions 
at the tip of the light-cone or at the intersection surface. In [11, Section 7.1] a scheme is 
presented where gi2\-yy is prescribed, together with wave-map conditions. It is obvious 
that there exist further schemes which are mixtures of the above and which might be 
more appropriate for some specific physical situations, or for matter fields with exotic 
coupling to gravity. 
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Rcndall's analysis, or that in [11], makes it clear that every vacuum characteristic 
data set as defined in Section 4.1 leads to a unique, up to isometry, associated space- 
time, either near the tip of the light-cone, or near the intersection surface S. Here 
uniqueness is understood locally, though again it is clear that unique maximal globally 
hyperbolic developments should exist in the current context. 

5. Solving the constraint equation 

There are several ways of solving (54). The aim of this section is to present those 
methods, in vacuum. One should keep in mind that some further specific hypotheses 
on the matter fields might have to be made in the schemes below for non- vacuum initial 
data: 

5.1. Solving for k 

For any g for which r has no zeros, (53)- (54) can be solved algebraically for k. This 
appears to be the most natural choice near the tip of a light-cone, where r is nowhere 
vanishing. 

5.2. T and [gjig] as free data 

Another way of solving (54) is to prescribe [gAs] ^^'^ the mean null extrinsic curvature 
T. Here one can simply choose r to be nowhere vanishing such that (54) is solvable for 
K. Regularity conditions on t in the light-cone-case are discussed in an appendix. 

5.3. K = 

Kendall's proposal is to reparamcterize the characteristic curves so that k — 0, (54) 
can then be rewritten as a linear equation for a conformal factor, say fi, such that 
gAB = ^"^lAB, where jab — [Vab] i^ freely prescribed; compare (7). 

5.4. K and [^ab] '*'' fr^^ data 

In some situations it might be convenient not to assume that k = 0, but retain a version 
of the conformal approach of Rendall. This requires only a few minor modifications in 
Section 2: ^suffices_to replace (6) by (19), (8) by (26), (11) by (33) and (10) by (32) 

with Tpjy, W'^' and W^ set to zero. As a matter of course the corresponding equations 
on N2 have to be adjusted analogously. 

By an appropriate choice of k, i.e. by choosing an adapted parameterization of the 
null rays, equation (19), which determines r, can sometimes be simplified; an example 
will be given in the next section. 

5.5. K = '''/(?"'■ ^ 1) 

An elegant approach is due to Hayward [27] who, in space dimension n = 3, proposes 
to use a paramcterisation where k — T/{n — 1). Then, in vacuum, (54) becomes a 
linear equation for r, in terms of the trace-free part of x which depends only upon the 
conformal class of g, 

diT + \a\^=0. (57) 

The solution r can then be used to determine a conformal factor O^ relating 'gAB with 
a freely prescribed representative jab of the conformal class, 

din - -^ (t - l-j^^'dijAB) = . 

n — 1 \ 2 / 
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In the case where data are given on a light-cone, one has to face the question of 
boundary conditions for (57), of the (necessary and/or sufficient) conditions on the 
data which will guarantee regularity at the vertex, and a possible relation between 
those. 

To address those questions, we start by comparing the k — t / [n — l)-gauge with 
the geometric k = Q gauge. Those quantities computed in the latter gauge will be 
labelled by ° in what follows. Both gauges are related by an angle-dependent rescaling 
of the coordinate r: Using the transformation law of the Christoffel symbols we find 

r/(n-l) = . = rh = |;f}, +|^|!; (58) 

or '^.^^ or Or'^ 

for r — r{f^x). That yields 

r{r)= f e^ Z"^' ^('■^^'^'■Mri = f e"^ ■'"■'■'■''' l'^('^=')l''^'"^'^'"Mri, (59) 

where we have suppressed any angle-dependence, and left unspecified any potential 
constants of integration. This defines the desired local diffeomorphism. 

As an example (and to obtain some intuition for this gauge scheme) consider the 
flat case where jcrp = and for which we can compute everything explicitly. There is 
no difficulty in determining the transformed data which satisfy jirp = 0. The general 
solution of (57) is 

T(r,x'^) = ro(x^) . 

Now we can explicitly compute (59), 

r(r)=A(i)+A(')e^'^ (60) 

for some integration functions A^^' , with A^'^' and tq nowhere vanishing since we seek 
a map r t-^ r which is a diffeomorphism on each generator. Then 

dr\ / / ,. N N n — 1 



Hr)=^-) r^ir)) = —^. (61) 

We choose, as usual, the affine parameter r in such a way that {r = 0} represents 
the vertex and such that f = ^^^. This leads to A^^^ = 0. Consequently, we either 
have to place the vertex at r = — oo and choose a positive tq or at r = -|-cx) with 
a negative tq (w.l.o.g. we shall prefer the first alternative). We conclude that in the 
K ~ T/{n — l)-gauge we need to prescribe initial data for aU r G M. 

The regularity condition for r translated into the k = t / {n — l)-gaugc docs not 
lead to any boundary conditions for r, except for the requirement of constant sign. It 
determines instead the position of the vertex, which in the new coordinates is located 
at infinity. 

Let us come back to the general case, which we tackle from the other side, namely 
by starting in the k ~ 0-gauge. We use the identity §p^ = — (§^) §^5^ to rewrite (58): 

dr T T f dr\ d^r 



dr n — \ 71 — 1 \dr J dr'^ 

d^r T dr 



dr^ n ^ I dr 

(r) = r e-^ /"' ^(^^^'^^^ dfi . (62) 



r 



This provides the inverse coordinate transformation. 

Now, for a smooth metric, in adapted null coordinates which are constructed 
starting from normal coordinates the generators are afiinely parameterized and we have 

T = ^^ + Oir) . (63) 
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But this behaviour remains unchanged under all reparameterisations of the generators 
which preserve the afBnc parameterisation as well as the position of the vertex. It 
follows that (63) holds for all smooth metrics in the gauge k = 0. 
From (62)-(63) we obtain 

r{f-) = A^i) + A^^'> logr + 0{r^) , A^^'> ^ Va;^ . 

If we start in the k — 0-gauge, with the vertex at f = 0, and transform into the 
K = T/{n — l)-gauge then, similarly to Minkowski space-time, the vertex is shifted to, 
w.l.o.g., r ~ — oo. Thus, space-time regularity forces the vertex to be located at infinity 
in the k = T/{n — l)-gauge. 

5. 6. The shear as free data 

Following Christodoulou [13], we let the second fundamental form x of a null 
hypersurface ,jV with null tangent i be defined as 

x{X,Y)^giVxe,Y), (64) 

where X,Y G T.A' . Choosing I to be dr we then have, using [11 , Appendix A], 



-M _ -fC 



XAB = gi^Adr, Ob) = 9^,B^Ar = 9cB^Ar + 9uB^ A 
_ 1 

~ 2 



7;drgAB > (65a) 



Xrr = g{Vrdr,dr) = , (656) 

XAr == g{S/Adr,dr) = 9 ^r^Ar = ffurCr = ■ (65c) 

Let a be the trace-free part of x on the level sets of r: 

t Qj^ 

CTAB := XAB -9 Xcd9ab ; 

71—1 

a is often called the shear tensor of -yV. It has been proposed (cf., e.g., [13]) to consider 
a as the free gravitational data at o/K. There is an apparent problem with this proposal, 
because to define a trace-free tensor one needs a conformal metric; but if a conformal 
class [g{r)] is given on ./K, there does not seem to be any need to supplement this class 
with a. This issue can be taken care of by working in a frame formalism, as follows: 

Let c/K be a n-dimensional manifold threaded by a family of curves, which we call 
characteristic curves, or generators. We assume moreover that each curve is equipped 
with a connection: if, in local coordinates, 9,. is tangent to the curves, then we let k 
denote the corresponding connection coefhcient, as in (51). 

Suppose, for the moment, that ^ is a characteristic hypersurface embedded as 
the submanifold {u = 0} in a space-time ^. Choose some local coordinates so that 
dr is tangent to the characteristic curves of ./K. Let Ca denote a basis of T^ along 
c/K such that 

VrCa = . (66) 

Let S denote an (n — l)-dimensional submanifold of ,JY (possibly, but not 
necessarily, its boundary) which intersects the generators transversally. We will further 
require on S that eo is null, and that for a ~ 2, . . . , n the family of vectors Ca is 
orthonormal. These properties will then hold along all those generators that meet S. 

Let x^ be local coordinates on 5, we propagate those along the generators of ,jy 
to a neighbourhood '^ C ^ of 5* by requiring S£q^x^ — 0. 

We choose ei ^ dr at S\ (51) implies then that this will hold throughout '^: 

ei = e/ar on ^. (67) 

We choose the Ca's, a = 2, ... n, to be tangent to S; since T,yV coincides with e-^, 
the Ca's, a = 2, . . . n will remain tangent to ^/K: 

Co = ^a dr + Ca^ds On 'W , a — 2, . . . ,n. (68) 
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On S wc choose the vector cq to be nuh, orthogonal to S, with 

5(eo,ei) = l. (69) 

Let {0°}a=o,i n, be a space-time coframe, defined on "^ C ^, dual to the frame 

{ea}a=o,i,...,n- From what has been said we have 

g = f (8)9^+9^ (E)e° + 9 ^ (8)9^ + ... + 0" 0'\ . (70) 

By construction, the one-forms 6*" are covariantly constant along the generators of ^: 

V^ei" = . (71) 

Here V is understood as the covariant-derivative operator acting on one-forms. 
Again by construction 0" annihilates T^, thus 0° ~ du along c/K: 

61" = 9°udu on '^. (72) 

We further note that 

r(a^) = Oon'^ fora==2,...,n. (73) 

To see this, recall that dr is orthogonal to Oa, hence 

N 
Q = g{dr,dA) = Tlab9''{dr)9\dB) ^Y.^'\dr)9%dB) • 

The result follows now from the fact that the (n — 1) x (n — 1) matrix {9°'{dB))a>2 is 
non-degenerate. 

We would like to calculate g := 'gAudx^dx^ using the covectors 9°^. For this, note 
that (70) and (72) imply 

N 

lAB^Y.^''{dA)e'^{dB). (74) 

a=2 

Thus, to determine 'g it sufhces to know the components {9°- b '■— 9°-{dB))a>2- Now, 
using (73) together with [11. Appendix A] we have for a > 2 

= V,J"s = drO'^B - T\Be\ = a,.r B - \g^^drgcB9''A . (75) 
There holds thus the following evolution equation for {9"^^ B)a>2'- 

dr^B - 5^^XCi3^"A = , (76) 

where 'g^'~' denotes the matrix inverse to X]a=2 9°" a9°' b- 
Let, as before, i = dr and define 

-B^. := V^4 . (77) 

Let Bab denote the frame-components of B: 

Bab = ea^efa'B^, ^=^ B,,, = 9-'{d^,)9\d,)Bab ■ (78) 

It follows from the definition (64) that x encodes the information on components of B 
in directions tangent to ^: 

XAB = Bab , Xrr = Brr = , XAr = Bat = . (79) 

The key distinction between B and x is, that B has components along directions 
transverse to ^, while x hasn't. Also note that for a,b> 1 the frame components Bab 
only involve the coordinate components of B^^ tangential to ./K, so the frame formula 

Xab = Bab , a,b>l, 
is geometrically sensible. 
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The last two equations in (79) give 

n n 

XAC = Y. ^''AO'cBab = Y. ^"AO^CXab , (80) 

a,b=2 a,b=2 

which allows us to rewrite (76) as 

1 = 0. (81) 



(82) 



(83) 

This equation leads naturally to the following pieture, assuming for simplicity 
vacuum Einstein equations. Consider, first, two null transversely intersecting 
hypersurfaces ^ and ^, with ^ n ^ = S. For a,b > 2 let 77°^ be one when 
a and b coincide, and zero otherwise. In addition to k, the gravitational data on 
^ = ^i U c/f2 can be encoded in a field of symmetric ?7-trace-free (n — 1) x (n — 1) 
matrices (Tati, a, & = 2, . . . , n. 

CTah = CTba , ?y" fab = . 

Let r be a solution of the equation 

(dr - h)t + + \a\l = , where \a\l := r^v'^'^abacd ■ (84) 

n — 1 ' ' 

There remains the freedom to prescribe r = g xab = X)" h=2 V'^^Xab on S* (one such 
function for each surface Ni and A^2)- Define 

T 
Xab = (Tab H -?7a6 ■ (85) 

n — I 
Solving (83) for 9°'b along the generators of -yPi and ^, we can calculate g^^ on ^ 
from (74), as long as the determinant of the matrix (6'°s)a>2 does not vanish (which 
will be the case in a neighbourhood of S). Here one has the freedom of prescribing 
0°'B on S' for a > 2. The characteristic constraint equation R^yl^^f = holds by 
construction. Indeed, the relation aab = Ca £^b oab, (iib > 2, can be justified in an 
analogous manner as equation (80). That gives, using (82), with a, b,c,d> 2, 

|a|2 = r^^'-rf'^aab'Jcd = 5^^e%e^c5''''^'s^''i3ea^e6^'T£f e.^ed^acH 
= g^^r'''TABacD^\a\\ 

and the assertion follows immediately. 

We can now apply any of the methods described previously (e.g., Kendall's original 
method if k = 0) to obtain a solution of the characteristic Cauchy problem to the future 
of^.ft 

One should keep in mind the following: prescribing the data Ca Is, or equivalently 
S°'b\s, determines the metric gAB\s on 5*. There is a supplementary freedom of making 
an 0{n — l)-rotation of the frame; 

ea%ix^) ^ io\{x^)et%{x^) , 

{t Note that prescribing a^b and k is equivalent to prescribe Xab as primary data. If we assume for 
simplicity that the r;-trace of Xab is nowhere- vanishing, one can then determine k from (84), and 
continue as described in the paragraph following (85). One could also consider this procedure in an 
adapted frame: From XAB on^ determines successively Qj^g, r and ft. However, since the 'gj^^-txzce 
of X is then not known a priori, it is not clear how to satisfy the constraint (84). 
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where the a;''a(x^)'s are 0{n — l)-matriccs. Any such rotation needs to be reflected in 
the (Jab's: 

aab{r,x'^) >^ uj''a{x'^)uj''b{x'^)crcd{r,x^) . 

So, in this construction aab undergoes gauge-transformations which arc constant along 
the generators, and are thus non-local in this sense. 

In the case of a light-cone, the above construction can be implemented by first 
choosing an orthonormal coframe 0° = (f)°' aAx^, a > 2, for the unit round metric 
SAsdx^dx^ on S"''^. The solutions 6''^ := re/)'' = rcjy'^Adx'^, a > 2, of (83) are then 
chosen as the unique solutions asymptotic to rcfia- It would be of interest to settle the 
question, ignored here, of sufficient and necessary conditions on Uab so that the resulting 
initial data on the light-cone can be realized by restricting a smooth space-time metric 
to the light-cone. 

5. 7. Friedrich 's free data 

In [20, 22], Friedrich proposes alternative initial data on ^/K, based on the identity ff 

dr'9AB - l^dr'gAB " -jf^" dr'g AC^rg BD = -^RatBt ] (86) 

equivalently 

QtBab — hBab — g BacBbd — —RatBt ■ (87) 

Equation (86) shows that the component RatBt of the Riemann tensor can be 
calculated in terms of n and the field 'gAB ■ 

Alternatively, given the fields CArBr, P = T^r and k, together with suitable 
boundary conditions, one can solve (86) to determine '^ab- ^Oj i^ space-time dimension 
four, Friedrich [20] proposes to use frame components of CatBt as the free data on ^ . 
There is, however, a problem, in that C ArBr is traceless 

= 3 (-^ firur = g ^ ArCr ■ 

So this condition has to be built-in into the formalism. But, as in the previous section, 
the tracelessness condition does not seemingly make sense unless the inverse metric 
g^^ , or at least its conformal class, are known. 

This issue can again be taken care of by a frame formalism, whatever the 
dimension, as follows: Let the orthonormal frame Ca, a = 0, . . . , n, and its dual coframe 
9°" be as in the last section. The property that the frame is parallel along the generators 
implies 

drea"" = -r^c^a^' = -g^'^BAcea^' for a > 2 . (88) 

We then have, for a,b> 2, 

drBab ^dr{ea^eb^B^^) 

= dr{ea^)eb^BAC + ea^dr{eb^)BAC + ea^Ck^drBAC 
-AD n E C u A-CD n E r, 

^ - g BoESa efc Bag - Sa g Boseb Bag 

+ Ca^Cb'^ {kBaC + V^^BaeBcD - RAtCt) 
= - ea'^Cb^-fBAGBoE + I^Bab - Ca^b'^RArBr ■ (89) 

Using 

n 

c,d=2 

ft This equation reduces to the usual Riccati equation (cf., e.g., [2()]) satisfied by the null extrinsic 
curvature tensor when k = 0. We are grateful to Jose-Maria Martin-Garcia for providing the general 
version of that equation. 
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we conclude that 

n 

drBab = — 2_^ ^"^ BacBdb + t^Bab ~ Gq Cfc RatBt ■ 
c,d=2 

From the definition of the Wcyl tensor in dimension n + 1, 



'^^uap • ^^ucrp 



n 



-i \9t^(y^i'p Qpp^ffj Qffj^pp I fjup^pa ) 



we find 



+ ^/ _ ^. R{9p<y9up - 9pp9va) , (90) 



CatBt — RatBt 79ABRr 

71 — 1 



For a, 6 > 2 let 



''Pab '■= Gq 66 CArBr 

represent the components of CatBt in the current frame. Then ipab is symmetric, with 
vanishing 77-trace. We finally obtain the following equation for Bab = Xab, a,b > 2: 

" 1 — 

(dr - K)Xab = - V" ■n'"^XacXdb " i>ab rVabTrr ■ (91) 

•^ — ' n — I 

c,d=2 

This equation shows that {Kjipab) can be used as the free data for the gravitational 
field: Indeed, given n, ipab and the component T^r of the energy-momentum tensor, we 
can integrate (91) to obtain Xab- Note that by taking the 77-tracc of (91) one recovers 
the constraint (84) (here with T„- possibly non- vanishing) . 

In the case of two transverse hypersurfaces the integration leaves the freedom of 
prescribing two tensors Xab on S*, one corresponding to A^i and another for A^2- Then 
one proceeds as in the previous section to construct the remaining data on the initial 
surfaces, keeping in mind the further freedom to choose 9°^B\s^ a > 2, on S. 

On a light-cone, (91) should be integrated with vanishing data at the vertex. 
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Appendix A. The expansion r and the location of the vertex 

The Hay ward gauge-condition of Section 5.5 has led us to the interesting conclusion, 
that in some gauge choices the vertex of the light-cone will be located at infinity. This 
raises the question: Under which conditions is the hypothesis, that the vertex is located 
at r = 0, consistent with natural boundary conditions at the tip of the light-cone? 

We start with the derivation of a necessary condition which needs to be imposed 
on the behaviour of the initial data in the k = t /(n— l)-gauge near the vertex in order 
to be compatible with regularity. It is known [LI] that in a k = 0-gauge arising from 
normal coordinates the initial data have to be of the form 

lAB = r'^SAB + hAB , where Hab = Oi{f^) . (A.l) 

Recall that we decorate with a circle the quantities corresponding to the gauge k = Q. 
We want to work out how such data look like in the n = T/{n — l)-gauge. The 
coordinate transformation (62), which defines a local diffeomorphism as long as A^'^'> 
does not change sign, reads 

r{f-) = A(i' + A(2) log r + fh A(2) ^ V x-^ 

= logr + A, (A.2) 
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where fh = Oi{r^) is determined by Hab and where A^'^' are integration functions. Here 
we have set A^^^ = and A^^^ = 1, so that the r-coordinate in the n = T/(n— l)-gauge 
is completely fixed, once f has been chosen. 

From (A. 2) we extract the behaviour of the inverse transformation, 

f{r)=e' + gh, gn = Oi{e^'') 
(where the symbol O in connection with the r-coordinate refers to the limit r — )• — oo). 
Now we can compute the overall form of the initial data, 

lAB{r) = lAB{r{r)) = e^^'sAB + ^As , where fc^s = Oi{e'^'') (A. 3) 
This implies 

W? = 4 (d^l^'d^lAB + tl^^^riAEl) = o(e-) . (A.4) 



4 \ n - 1 

Note that, in contrast to the k — 0-gauge, r remains bounded at the vertex for 
regular light-cone data of the form (A. 3), due to (57) and (A.4). 

Next, let us show that a bounded r can only be compatible with regularity when 
the vertex is located at infinity. For definiteness, we consider initial data g^^ with 
nowhere vanishing r = ^S ^iSabj within the scheme of Section 3. Then k is 
computed algebraically via (21) (and depends on the initial data). Note that at this 
stage r is a known function of r which can be regarded as "gauge part" of the initial 
data. 

By calculations similar to those in (58)- (59) we can then obtain the coordinate f 
relevant to the k — 0-gauge: 

r{r) = I e^"' <^--)'^^^dri , (A.5) 

and transform all the fields to this gauge. 

We have the identity 

or Ti — 1 

T(r) = -r-T(f(r)) , where f = — ; h 0{r) , (A. 6) 

or r 

since we assume regular light-cone data. We consider the maximal range of f, near 
r = 0, where f is positive. It follows from (A.5) that r H> r{r) is monotone 
there. Let us assume that this function is strictly increasing (the decreasing case 
is handled in a similar way), and let (i?i,i?2) denote the corresponding range of r, 
with — oo < Ri < R2 < oo. Then r is positive on (i?i,i?2) by (A. 6). If we choose 
i?i < ro < R2 such that r^^(ro) = fo > 0, from the last equation we find, for some 
(a; ^-dependent) constant A, 

rdf ^ / — = h 0{f) 1 dr = ^ + (n - 1) logr + 0{r'^) . 

Ira J ro \ ^ / 

The right-hand side diverges to minus infinity at the vertex f = 0, which is mapped to 
Ri. This gives 

rra 

rdf = +00 . (A. 7) 

We conclude that any gauge in which r is bounded will force the vertex to lie at infinity, 
for initial data which can be realized by a smooth space-time metric. 

As another application of (A. 7) we reconsider the k = T/(n — l)-gauge. Let us 
denote by ro(x ) the integration function which arises in the associated constraint 
equation diT+ jtrp = 0. The exponential decay of |crp at a regular vertex, cf. equation 
(A.4), is compatible with (A. 7) only if tq is bounded away from zero. 

Finally, consider initial data on (0,oo) with 

9ab =-r^SAB + Oi(r'*) , 
as in Section 3.2. The function r satisfies then 

T = ^^^ + 0{r) , (A.8) 

r 

which is, not unexpectedly, fully compatible with (A. 7). 
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